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PREFACE 

This  Memorandum  is  part  of  Rand's  study  for  the  Advanced  Research 
Projects  Agency  of  those  phenomena  which  affect  the  performance  of 
optical  reconnaissance  and  guidance  equipment.  The  objective  of  these 
studies  is  to  provide  sufficient  understanding  to  permit  the  systems 
analyst  to  compute  performance  estimates  under  various  operational 
conditions . 

In  optical  communications  and  related  devices,  the  random  variations 
In  the  received  signal  due  to  atmospheric  turbulence  can  represent  a 
severe  limitation  to  system  performance.  Theoretical  studies  of  these 
fluctuations  have  been  based  on  solutions  to  the  wave  propagation 
equation  which  are  correct  only  to  first  order  In  the  refractive  Index 
deviation.  It  Is  shown  in  this  Memorandum  that  a  solution  correct  to 
all  orders  Is  readily  obtained  by  direct  analogy  with  a  method  commonly 
adopted  to  solve  the  Schrodinger  equation  for  high  energy  potential 
scattering.  Study  of  the  more  highly  developed  scattering  theory  provides 
useful  insights  Into  the  nature  and  limitations  of  the  various  approxi¬ 
mate  solutions  of  propagation  theory;  in  particular,  the  Memorandum 
points  out  that  the  standard  approximation  adopted  in  propagation 
studies  is  based  on  an  approach  known  to  be  invalid  in  scattering 
theory.  The  results  should  be  of  interest  to  those  concerned  with 
the  use  of  lasers  in  the  atmosphere. 
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SUMMARY 


By  direct  analogy  with  a  method  developed  by  Schiff  to  solve  the 
Schrddinger  equation  for  high  energy  potential  scattering,  it  is  pos¬ 
sible  to  solve  the  equation  for  wave  propagation  in  a  turbulent  medium 
in  a  manner  which  explicitly  demonstrates  that  the  solution  so  obtained 
is  correct  to  all  orders  in  the  refractive  index  deviation  and  to  low¬ 
est  order  in  the  stationary  phase  approximation.  Although  the  solution 
is  readily  extended  to  next  order  in  atationary  phase,  such  an  extension 
is  recognized  in  scattering  theory  as  unwarranted  since  it  neglects 
terms  of  the  sa  <«*  order  from  outside  the  region  of  stationary  phase. 

The  conventional  "Bom"  and  "RytovM  solutions  widely  adopted  in  propa¬ 
gation  theory  are  of  questionable  validity  since  they  represent  approx¬ 
imations  (first  order  in  refractive  index  deviation)  to  the  extended 


solution 
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THE  MOLIERE  APPROXIMATION  FOR  WAVE  PROPAGATION  IN  TURBULENT  MEDIA 

It  is  the  purpose  of  this  note  to  point  out  the  applicability  of 
the  Moliere  approximation,  used  in  high  energy  potential  scattering 
theory,  to  treatments  of  wave  propagation  in  random  or  turbulent  media. 
In  its  usual  formulation,  the  Moliere  approximation  to  the  solution  of 
the  SchrBdinger  equation 


is 


(1-4) 


[V2  +  k2  -  U(r)]  t(r)  -  0 


(1) 


2 

♦  ~  tM  “  exp  [* **kz  -  2k  [  U(x»y»z')  dz']  (2) 

•■00 
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where  e  represents  the  Incident  wave  function.  The  integration 

* 

indicated  in  Eq.  (2)  is  along  a  straight  line  parallel  to  the  z-axis. 
With  the  condition  ka  »  1,  where  the  length  a  represents  the  character¬ 
istic  distance  over  which  U(r)  changes  by  a  significant  fraction  of 
itself,  the  approximation  applies  throughout  and  near  the  region  of 
the  scattering  potential,  and  may  be  used  in  the  standard  integral  ex¬ 
pression  to  determine  the  scattering  amplitude.  Although  limited  to 

**  , 

determination  of  small  angle  scattering,  it  may  be  used  for  strong 
potentials  where  the  Born  approximation  is  of  little  value. 

The  wave  equation  for  a  scalar  wave  propagating  in  a  medium  whose 
index  of  refraction  is  n(r) , 

*The  WKB  approximation  is  f  ~  tM  ■  exp  (i  f Vk2  -  U  ds),  with  the 
Integration  along  the  classical  trajectory. 

**SchiffO)  has  extended  the  Moliere  approximation  for  use  in 
large  angle  scattering. 


i 


2 


(V2  +  n2k2)  * **  -  0  (3) 

may  be  made  formally  equivalent  to  Eq.  (1)  with  the  substitution 

U(r)  ■  -  ek2,  where  e(r)  ■  n2  -  1. 

With  the  assumption  that  n  ■  1  when  z  <  0,  the  Moliere  approxi- 
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mation  for  the  wave  in  the  "near  zone"  (0  <  z  ~  ka  )  may  be  expressed 
in  the  form 


iS 


♦m  ■  *„• 


ikz  * 

where  *  ■  e  is  the  incident  wave,  and 

To  ’ 


(4) 


S 

o 


«(x,y,*')  dz' 


(5) 


Expansion  of  Eq.  (A)  gives  the  Born  approximation, 

*B  •  *o(l  +  1S„>  (6) 

In  treatments  of  propagation  in  random  media,  the  index  of  re- 

irk 

fraction  is  usually  expressed  in  the  form  n(r)  ■  1  +  npr) ,  hence 
2 

e  ■  2nj  +  n^,  and  Sq  may  be  expressed  in  the  form 

SQ  «  k  n1(x,y,z/)  dz'  +  j  J  n2(x,y,z')  dz '  (7) 


l 


*In  the  WKB  approximation,  *  ■  exp  ( tk  J*  nds),  where  the  inte 

gration  is  along  the  ray  path.  w 

** 

It  will  be  assumed  that  the  mean  value  of  n^  satisfies 
(n^)"  0,  and  that  (np  is  constant. 
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If,  as  is  often  assumed,  |n^|  «  1,  the  Born  approximation  may  be 
expressed  as 

♦B  "  ♦o(l  +  iSo>  (8) 


where 


The  corresponding  "Rytov  approximation" 


(9) 


(10) 


is  somewhat  ambiguous  in  that  its  exponential  form  implies  a  higher 

order  (in  n^)  approximation. 

(3) 

Schiff  has  shown  explicitly  how  the  infinite  Born  expansion 
may  be  summed,  with  each  term  in  the  expansion  evaluated  to  lowest 
order  in  the  stationary  phase  approximation,  to  give  the  Mollere 
approximation.  Thus,  Eq.  (4)  gives  the  wave  in  the  zeroth “order 
stationary  phase  approximation  correct  to  all  orders  in  n^;  the  Born 
and  Rytov  approximations,  Eqs.  (8)  and  (10),  are  zeroth-order  stationary 
phase  approximations  correct  to  first  order  in  n^.  Schiff  points  out 
that  the  next  order  stationary  phase  approximation,  given  by 


♦o« 


IS  +S. 
o  1 


(11) 


with 
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1  cz 

S1  “  -  4  Jo  (z  -  z')  h(x,y,z')  dz '  (12) 

2 

where  h  ■  7  e,  Is  Inconsistent  in  that  it  neglects  terms  from  outside 

the  stationary  phase  region  which,  in  general,  are  of  the  same  order. 

The  stationary  phase  expansion  is  an  expansion  in  powers  of  (ka) 

(2) 

Gol'dman  and  Migdalv  '  have  emphasized  that  such  "classical"  ex¬ 
pansions  are  asymptotic  and  can  never  give  diffraction  effects  cor¬ 
rectly. 

Similarly,  the  first-order  (in  n^)  approximation  to  Eq.  (11), 
in  either  the  Born  or  Rytov  form, 

“  V1  +  iSo  +  Sl)  (13) 


or 


iS'+S' 
o  1 


(14) 


where 


S'  -  -  £  J  (z  -  z')  h'(x,y,z ')  dz  *  (15) 

2 

with  h '  ■  27  n^,  is  also  inconsistent. 

Equations  (13)  and  (14)  have  been  widely  adopted  in  the  propaga¬ 
tion  literature  as  approximations  to  the  wave  in  the  near  zone.  Their 
inconsistency  is  not  manifest  in  the  conventional  development 
where  the  wave  equation  is  expanded  in  powers  of  n^  with  Sq  and 
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treated  as  quantities  o£  the  same  order.  Only  after  the  first  order, 
or  Rytov,  equation  is  reformulated  as  an  integral  equation  is  the 
stationary  phase  approximation  introduced,  essentially  by  restriction 
of  the  region  of  integration  to  the  "Fresnel  cone."  This  cone  has 
its  vertex  at  the  field  point  and  opens  in  the  negative  z-direction 
with  angular  spread  ~  l/ka;  the  same  angle  characterizes  the  "dif¬ 
fraction  cone,"  which  has  its  vertex  at  the  scattering  inhomogeneity, 
or  "turbule,"  and  opens  in  the  positive  z-direction.  The  Fresnel 
approximation  is  introduced  as  a  consequence  of  the  well-established 
"forward  scattering"  approximation  (the  assumption  that  essentially 
all  the  diffracted  radiation  is  confined  to  the  diffraction  cone) 
and  is  conjectured  to  be  of  general  validity  provided  ka  »  1. 

This  view  contrasts  strongly  with  that  adopted  in  potential  scattering 
theory,  where  the  use  of  the  Fresnel  approximation  is  restricted  to 

the  determination  of  S  . 

o 

It  is  of  interest  in  this  regard  to  consider  a  simple  example, 
the  scattering  of  a  plane  wave  by  a  single  spherical  turbule  with 
radius  a  and  constant  index  deviation  n^  with  the  condition  |n^ka|  «  1. 
As  is  well  known,  the  scattered  field  may,  with  reasonable  unambiguity, 
be  separated  into  a  widespread  (i.e.,  throughout  distances  of  order  a 
in  the  x,y-plane),  weak  coherent  "refractive"  field  and  a  narrow, 
intense  diffraction  beam.  If  a  field  point  very  near  (but  not  neces¬ 
sarily  on  the  axis  of)  the  turbule  is  considered,  that  portion  of  the 
turbule  within  the  Fresnel  cone  gives  the  refraction  field  with  excel¬ 
lent  accuracy,  but  gives  no  diffraction  field.  At  a  larger  distance, 

*  (3  7) 

More  precisely,  the  Fresnel  paraboloid  of  revolution  ’  (the 

region  of  stationary  phase) . 
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the  Fresnel  cone  Intersects  a  greater  section  of  the  turbule,  and 

implies  the  existence  of  a  weak,  wide-angle  diffraction  beam,  in 

addition  to  the  refraction  field.  As  the  distance  is  increased,  the 

diffraction  pattern  as  determined  by  the  material  in  the  pertinent 

Fresnel  cone  increases  in  intensity  (and  total  scattered  energy)  and 

decreases  in  angular  opening;  the  refraction  fields  remain  accurately 

2 

determined.  For  distances  exceeding  ka  ,  the  Fresnel  cones  associated 
with  points  in  the  scattered  field  intercept  essentially  the  entire 
turbule,  and  the  approximation  is  obviously  valid. 

The  errors  associated  with  this  initial  elimination  and  gradual 
mixing-in  of  the  diffraction  field  depend  on  the  quantities  being 
determined.  For  example,  the  refractive  and  diffractive  phase  shifts 
are  of  the  same  order,  and  the  narrowness  of  the  diffraction  beam 
insures  that  the  near-zone  transverse  mean-square  phase  deviation 
associated  with  a  collection  of  independent  spheres  with  randomly 
distributed  values  of  n^  will  have  negligible  error  in  the  Fresnel 
approximation.  In  contrast,  the  amplitude  perturbation  in  the 
diffraction  field  is  considerably  larger  than  that  in  the  refraction 
field,  and  it  is  easy  to  show  that,  despite  the  narrowness  of  the  dif¬ 
fraction  field,  the  amplitude  fluctuation  associated  with  scattering  by 
a  densely  packed  collection  of  random  spheres  is  everywhere  dominated 
by  the  diffractive  effects.  The  example,  while  highly  artificial,  is 
in  agreement  with  the  general  conclusions  of  potential  scattering 
where  the  Fresnel  approximation  is  used  to  determine  Sq  (a  phase 
perturbation)  and  where  (the  Fresnel  amplitude  perturbation)  is 
regarded  as  inconsistent. 
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Nevertheless,  the  admonition  of  Schiff  and  others  against  the 
extended  use  of  the  Fresnel  approximation  is  of  necessity  somewhat 
imprecise  in  that  it  is  based  on  the  observation  that  the  contri¬ 
butions  from  outside  the  Fresnel  cone  are  likely  to  be  of  the  same 
order  as  S  ;  it  remains  possible  that,  under  certain  circumstances, 
the  outside  contributions  may  be  considerably  smaller  than  those  from 
inside.  Although  use  >_f  the  general  wave  equation  to  estimate  con¬ 
ditions  for  which  the  extended  Fresnel  approximation  may  be  valid 

presents  a  formidable  task,  it  is  possible  to  show'  '  that  when 
2  A  5 

z/(n^)  k  a  >1,  the  conventional  near-zone  Fresnel  solution  (Eq .  (14)) 

of  the  Rytov  equation  agrees  with  the  assumptions  underlying  the 

derivation  of  the  Rytov  equation  as  an  approximation  to  the  wave 

equation.  While  this  "consistency"  condition,  based  on  a  solution 

of  the  approximate  equrtion,  is  neither  necessary  nor  sufficient  for 

the  validity  of  the  extended  Fresnel  approximation,  such  conditions 

are  commonly  adopted  when  more  rigorous  conditions  (those  based  on 

solutions  of  the  general  wave  equation)  are  unavailable. 

A  more  familiar  consistency  condition  on  the  Rytov  approximation 

is  the  Born  condition  |Sq|  <  1,  which  implies  (S^)  <  1;  and  hence  (note 

the  Fresnel  approximation  is  valid  for  the  determination  of  S  )  the 

o 

2  2  * 

"upper  limit"  condition  (n^)  k  az  <  1.  Thus,  values  of  z  for  which 

the  conventional  Fresnel-Rytov  solutions  are  consistent  may  be 

2  3 

determined  from  the  parameter  -y  ■  (n^)(ka)  :  When  y  <  1  the  solutions 

“ft 

This  condition  is  neither  necessary  nor  sufficient  for  the 
validity  of  the  Rytov  equation  and  has  been  derived  many  times  with 
various  approaches,  some  of  which  are  discussed  in  Ref.  8. 
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are  consistent  (although  not  necessarily  valid)  in  the  interval 
2 

y  £  z/ka  <  y  ;  when  y  >  1  the  solutions  are  inconsistent  (although 
not  necessarily  invalid)  for  all  z.  Reasonable  estimates  for  the 
various  parameters  indicate  y  »  1  for  optical  propagation  in  the 
atmosphere,  hence,  in  this  application,  the  Fresnel-Rytov  solutions 
should  be  used  with  caution. 


CONCLUSIONS 

Schiff's  approach  to  the  solution  of  the  Srhrodinger  equation 

(3) 

for  high  energy  potential  scattering  may  be  used  to  solve  the 

scalar  wave  equation  for  propagation  in  a  turbulent  medium  in  a 

manner  which  explicitly  demonstrates  that  the  Moliere  solution 

(Eq.  (4))  is  correct  to  all  orders  in  the  refractive  index  deviation. 

The  method  is  readily  extended  to  higher  order  in  stationary  phase, 

(2  3) 

but  the  extended  solution  (Eq.  (11))  is  invalid  ’  since  it 
neglects  terms  of  the  same  order  from  outside  the  stationary  phase 
region.  The  conventional  near -zone  Born  and  Rytov  solutions  are 
approximations  (of  order  n^)  to  the  extended  solution,  and  their 
use  is  unjustified. 


A  closely  related  inconsistency  is  that  which  arises  from  the 
assertion  (Refs.  5  and  6)  that,  for  z  ~  ka2 ,  the  Rytov  equation  may 
be  further  approximated  and  render  valid  the  equations  of  linear 
(in  n.)  geometric  optics.  This  view,  when  combined  with  the  widely 
adopted  upper  limit  given  above  (yz  <  ka2)  implies,  inconsistently, 
the  existence  when  y  >  1  of  an  interval  Y"*  <  z/ka 2  <  1  in  which  the 
Rytov  equation  is  Invalid  and  the  less  general  geometric  equation  is 
valid . 
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